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ABSTRACT 

Observations of the lio — loi 556.936 GHz transition of ortho-water with the 
Submillimeter Wave Astronomy Satellite (SWAS) have revealed the presence of 
widespread emission and absorption by water vapor around the strong submil- 
limeter continuum source Sagittarius B2. An incompletely-sampled spectral line 
map of a region of size 26 x 19 arcmin around Sgr B2 reveals three notewor- 
thy features. First, absorption by foreground water vapor is detectable at local 
standard-of-rest (LSR) velocities in the range -100 to km/s at almost every 
observed position. Second, spatially-extended emission by water is detectable at 
LSR velocities in the range 80 to 120 km/s at almost every observed position. 
This emission is attributable to the 180-pc molecular ring identified from previ- 
ous observations of CO. The typical peak antenna temperature of 0.075 K for 
this component implies a typical water abundance of 1.2 - 8 xlO~^ relative to 
H2. Third, strong absorption by water is observed within 5 arcmin of Sgr B2 
at LSR velocities in the range 60 to 82 km/s. An analysis of this absorption 
yields a Hg^O abundance ~ 2 — 4 x 10^^ relative to H2 if the absorbing water 
vapor is located within the core of Sgr B2 itself; or, alternatively, a water column 
density of ~ 2.5 — 4 x 10^^ cm^^ if the water absorption originates in the warm, 
foreground layer of gas proposed previously as the origin of ammonia absorption 
observed toward Sgr B2. 

Subject headings: ISM: abundances - ISM: molecules - ISM: clouds - molecular 
processes - radio lines: ISM - submillimeter 
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1. Introduction 

The Sagittarius B2 cloud represents the largest condensation of molecular gas and dust 
in the Galaxy. With a total gas mass of ~ 4 x IO^Mq (Lis & Goldsmith 1989) concentrated 
within a region of size ~ 30 pc, Sgr B2 exhibits the largest submillimeter continuum flux of 
any astronomical source outside the solar system. 

The molecular inventory of the gas in Sgr B2 has been extensively studied by means of 
ground-based radio observations (e.g. Cummins, Thaddeus & Linke 1986; Turner 1989, 1991; 
Sutton et al. 1991; Nummelin et al. 1998, 2000). Such observations have been supplemented 
by far-infrared and submillimeter observations carried out from airplane altitude or from 
space. In particular, observations carried out over the past decade with the Kuiper Airborne 
Observatory (KAO; see Zmuidzinas et al. 1995), the Infrared Space Observatory (ISO; see 
Cernicharo et al. 1997) and the Submillimeter Wave Astronomy Satellite (SWAS; see Neufeld 
et al. 2000) have facilitated the observation of water vapor in Sgr B2. Water is of particular 
interest because of its central role in the chemistry of interstellar oxygen and its potential 
importance as a coolant of molecular clouds (see, for example, Neufeld, Lepp & Melnick 
1995; Melnick et al. 2000; Bergin et al. 2000). 

In a earlier study of Sgr B2, we used SWAS to obtain the first high-spectral-resolution 
[Av ~ Ikms^^) observations of Hg^O toward that source. In a previous paper (Neufeld et 
al. 2000; hereafter NOO), we reported the detection of water absorption over a wide range 
of LSR velocities that can be attributed to water vapor in foreground clouds that are not 
associated with Sgr B2. We have since extended these observations to map (incompletely) a 
region of size 26' x 19' around Sgr B2. In the present paper, we describe these new mapping 
observations and discuss the abundance of water vapor in the vicinity of Sgr B2 itself. 

2. Observations 

Our mapping observations of Sagittarius B2 were carried out during the period 1999 
March 20 through 2000 September 16. They yielded incomplete coverage of a region ex- 
tending 9.6' North, South and East and 16.0' West of Sagittarius B2. In Table 1 we list 
the on-source integration times for each of 19 positions observed, together with the offsets 
relative to the core of Sgr B2 at position a = 17^47" 19.^7, S = -28°23'07" (J2000). At 550 
GHz the SWAS beam size is 3.3' x 4.5' and the main beam efficiency is 0.9 (Melnick et al. 
2000). All the data were acquired in standard nodded observations (Melnick et al. 2000) and 
were reduced using the standard SWAS pipeline. The reference (nod) position was located 
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at a = 17^54"' 12.^2, 5 = -29°42'19" (J2000) in a region of devoid of strong CO emission^. 
The total on-source integration time for these mapping observations was 90.4 hours. 

Because SWAS employs double-sideband receivers, the ^^CO J = 5 — 4 transition at 
550.927 GHz and the H2O lio — loi transition at 556.936 GHz are observed simultaneously 
in the lower and upper sidebands of one receiver (Receiver 2). For the mapping observations, 
the SWAS local oscillator (LO) setting was chosen to separate the ^'^CO emission from Sgr 
B2 (at fLSR ~ 65kms^^) from any H2O emission or absorption at LSR velocities Vlsr ^ 0. 
For the (0,0) position, several additional LO settings had previously been used (see NOO) to 
allow the region —145 km s~^ < f lsr ^ km s~^ to be covered as well without contamination 
by ^'^CO emission from the lower sideband. 



3. Results 

In Figure 1, we show the nineteen H2O spectra obtained from our mapping observations 
described above. Each panel shows a spectrum covering the range —145 km s^^ < ^lsr ^ 
145kms~^. The vertical dotted line denotes ^lsr ~ 65kms^"'^, the systemic velocity of Sgr 
B2. The upper horizontal dotted lines indicate the continuum flux level, determined from 
the spectral regions vlsr = —145 to — 125kms~^ and 125 to 145kms~^ that are devoid of 
absorbing gas, and the lower horizontal dotted lines indicate the zero flux level, estimated 
by the method described in §4.1 below. On each panel a pair of numbers, X/Y, describe 
respectively the continuum antenna temperature (which corresponds to the separation be- 
tween the two dotted horizontal lines), and the antenna temperature (in K) at the top of 
the panel. The former is also given in Table 1. 

The feature at i'lsr ~ —90 to — 50kms~^ is a ■'^^CO emission feature detected in the 
lower sideband. In Figure 2, we show the Receiver 2 data reduced as appropriate for the 
^'^CO line in the lower sideband; a strong feature at ^lsr ~ 65kms"^ is observed at several 
positions close to Sgr B2. We note that the spectra in Figure 2 arc not perfect mirror images 
of those shown in Figure 1; in particular, the ^"^CO features are significantly narrower in 
the spectra that were reduced as appropriate for the lower sideband (Figure 2). The reason 
for this behavior is that the SWAS LO setting is fixed in the rest frame of the spacecraft, 
which has a motion relative to the local standard of rest (LSR) that varies quite rapidly. 
In particular, the spacecraft orbits the Earth every 97 minutes with an orbital speed of 
7kms~^ (and - of course - the Earth also orbits the Sun.) In the standard SWAS pipeline, 
the necessary Doppler corrections are applied to every 2 seconds' worth of data, and the 



At this reference position, the antenna temperature for ^^CO J=l-0 was less than 0.5 K at all velocities 
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sign of the correction is opposite for the lower and upper sidebands. Thus any spectral line 
originating in the lower sideband will be smeared out as it appears in Figure 1, while any 
line from the upper sideband will appear broadened in Figure 2. 

Although not the primary emphasis of this paper, for completeness we present additional 
spectra obtained for the ^Pi — ^ Pq fine structure line of atomic carbon. These spectra, shown 
in Figure 3, are obtained from Receiver 1 (upper sideband) on the SWAS instrument, which 
performs simultaneous observations of the CI ^Pi — "^Pq line in a linear polarization orthogonal 
to that observed with Receiver 2. There is no evidence for emission or absorption by the 
487.249 3i — 82 transition of O2, to which Receiver 1 is sensitive in its lower sideband. 

4. Discussion 

Even a cursory inspection of Figure 1 reveals three significant features in the water spec- 
tra: (1) extended absorption, observed at almost every one of the 19 positions, covering the 
range — lOOkms"^ ^ "^lsr ~ 0; (2) extended emission, again observed at almost every posi- 
tion, covering the range 80 ^ t'LSR ~ 120kms~^; (3) less-extended absorption, only observed 
unequivocally at the (0.0', 0.0'), (0.0, +3.2'), (0.0,-3.2'), (-3.2', 0.0) and (-3.2', -3.2') po- 
sitions, centered at vlsr ~ 65kms~^, the systemic velocity of the Sgr B2 cloud. 

4.1. Extended absorption by water vapor at vlsr ~ Oknis"^ 

Broad absorption by water vapor over the entire velocity range wlsr = ^100 to Okms^^ 
has been inferred previously from SWAS observations of H2^0 and Hg^O toward Sgr B2 
(NOO; see also Figure 7), the H2^0 spectrum indicating that the peak optical depth occurs 
s-t vi^SR ~ 0. The absorption in this velocity range can be attributed to material lying at a 
wide range of locations along the hne-of-sight between the Solar circle and Sgr B2. 

The mapping observations presented in Figure 1 show clear evidence for absorption 
at Vlsr ~ Okms~^ at 18 of the 19 positions observed, implying that both the continuum 
emission and intervening material have a large angular extent. The angular extent of the 
absorbing water vapor at vlsr ~ is considerably larger than that of the water vapor at 
fLSR ~ 65kms~^, the systemic velocity of Sgr B2. The large spatial extent of the water 
absorption at vlsr ~ Okms~^ is further demonstrated by the SWAS spectrum of water 
vapor toward Sgr A, shown in Figure 4, which clearly exhibits strong absorption in this 
velocity range. Located at position a = 17^' 45" 39.^^7, 6 = -29°00'30" (J2000), SgrA is 
at offset (-21.8', -37.4') relative to Sgr B2. Thus the SWAS observations of the Galactic 
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Center region show water absorption at vi,sr ~ km s~^ covering a region of longest angular 
dimension 50 arcmin, corresponding to a linear scale of 125 pc for an assumed distance to 
the Galactic Center of 8.5 Kpc. Of course, if the absorbing material lies significantly closer 
to us than the Galactic Center, then its linear extent could be smaller than 125 pc. 

The presence of strong absorption at vi,sr ~ allows us to estimate the continuum flux 
level, a quantity that is not reliably determined by the SWAS instrument. In particular, the 
zero flux level in Figure 1 and Table 1 is established on the assumption that the spectral 
region Wlsr = ^5 to 5kms^^ is subject to complete absorption by water vapor in foreground 
material along the linc-of-sight to Sgr B2. This assumption is motivated by the facts that 
(1) absorption at wlsr ~ Okms^^ is observed at almost every position and is well separated 
from contamination from ^^CO emission in the lower sideband; and (2) the H2^0 spectrum 
of the lio — loi hue at the (0.0', 0.0') position (see NOO; also Figure 7) implies an optical 
depth ~ 1 at vlsr ~ 0, corresponding to an H2^0 lio — loi optical depth of several hundred. 
If the I'LSR = ~5 to 5kms~^ spectral region is not subject to complete absorption, then the 
continuum antenna temperatures quoted are lower limits. 

The continuum antenna temperatures obtained by this method can be compared with 
continuum maps presented by Dowell et al. (1999), who observed Sgr B2 with the SHARC 
instrument at the Caltcch Submillimcter Observatory. In Table 1, we show the 350 fim 
continuum fluxes (in kJy per SWAS beam) - obtained by convolving the SHARC data with 
the SWAS beam profllc^ - along with the SWAS-determined 550 GHz ( 540 /im) continuum 
fluxes. The SHARC-mcasurcd 350 fim continuum fluxes should be regarded as lower limits 
because the relatively small offset of only 2 arcmin for the SHARC reference beam (Dowell, 
private communication) can lead to highly extended emission being "chopped out"^. Toward 
the central 9 positions [i.e. the (0,0) position and 8 positions closest to it], the inferred 540 
lim / 350 fjLm continuum flux ratios all he in the range 0.17 to 0.42. We have compared these 
values with those expected for optically-thin dust with an opacity that varies as A~^. The 
ratios observed toward the central nine positions are all consistent with dust temperatures 
of 15 K and above; the larger values of the 540 fim. / 350 /im flux ratio observed at greater 
distances from the cloud center suggest that the dust temperature is smaller in the outer 
parts of the Sgr B2 cloud. 



'"'approximated as a Gaussian circular beam of diameter 4 arcmin (FWHM) 

^ A further uncertainty in the comparison of the SWAS and SHARC continuum measurements results from 
the different assumptions used in calibrating these measurements. The SHARC data calibration effectively 
assumes that the beam is coupled to a point source, while the SWAS data is calibrated assuming that the 
source fills the main beam of the antenna. 
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4.2. Extended emission from water vapor in the 180-pc moleculcir ring 

(vlsr ~ 80 - 120kms-i) 

Extended water emission in the Wlsr ~ 80 — 120kms~^ velocity range is detected from 
18 of the 19 positions we observed in the vicinity of Sgr B2, with the integrated antenna 
temperatures given in Table 2. Emission in this velocity range is also evident in the spectrum 
of Sgr A (Figure 4), although it is blended with a strong emission feature centered at i'lsr ~ 
65kms~^. In Figure 5, we show the average spectrum for those 14 positions near Sgr B2 
that are devoid of apparent absorption at the Sgr B2 systemic velocity. The solid line is 
Gaussian fit to the average emission feature: its velocity centroid is 93kms~^, the line width 
is 43kms^^, the peak antenna temperature is 0.075 K, and the velocity- integrated antenna 
temperature is 3.4Kkms~^. 

This extended emission feature is kinematically similar to a highly extended emission 
component detected in ^'^CO maps of the Galactic Center region (Bally et al. 1987; Lis & 
Goldsmith 1989, hereafter LG89). It arises from the far side of a molecular ring of angular 
diameter ~ 2 degrees - corresponding to a projected radius of ~ 180 pc - the near side 
of which shows ^^CO emission and H2O absorption (NOO) at negative LSR velocities. The 
kinematic properties of this 180-pc molecular ring yield a parallelogram in the I — v plane; 
they have been attributed either to an explosion-driven expansion (Scoville 1972; Kaifu et 
al. 1972) or - more recently - to the gravitational effects of a Galactic bar potential (Binney 
et al. 1991). 

In Table 2, we also show the H2O lio — loi / ^^CO J = 1 — line ratio for each of the 
SWAS observed positions. The ^^CO J = 1 — spectra were obtained by convolving the map 
of Bally et al. (1987) with the SWAS beam profile, approximated as a Gaussian circular beam 
of diameter 4 arcmin (FWHM). For the 14 positions that are devoid of apparent absorption at 
the Sgr B2 systemic velocity, the line ratios apply to the velocity range vlsr — 80—120 km s~^ . 
For the other 5 observed positions, the line ratios apply only to emission redward of the 
average H2O line centroid, i.e. to the velocity range ^lsr = 93 — 120 km s^^. The reason for 
restricting the velocity range for these positions is to avoid any effects of the Sgr B2 water 
absorption line upon the line ratios derived for the extended emission feature. 

With the exception of a single position at which the H2O lio — loi line was not detected, 
the H2O lio ~ loi / ^^CO J = 1 — line ratio lies within a factor 2 of 0.065. In order to 
estimate the water abundances implied by the observed II2O lio — loi 

/ 13C0 J= 1-0 line 

ratios, we have made use of a statistical equilibrium calculation for the line strengths. We 
adopted the rate coefficients of Flower & Launay (1985) for coUisional excitation of CO by 
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H2, and those of Phillips, Maluendes & Green (1996) for that of H2O by H2/ Our excitation 
model indicates that the ^^CO J = 1 — line is optically thin and that the H2O lio — loi 
is effectively thin*^. Thus the inferred H2O to ^^CO abundance ratio is proportional to the 
H2O lio - loi / ^^CO J = 1 - line ratio. 

The inferred H2O to ^^CO abundance ratio also depends upon the conditions of temper- 
ature and density assumed for the emitting region. Based upon multitransition observations 
of CO using the AST/RO telescope, Kim et al. (2002) have argued that the physical condi- 
tions within a few hundred parsec of the Galactic Center are remarkably constant - except 
within the Sgr A and Sgr B regions themselves - with the gas temperature everywhere in 
the range 30 to 45 K and the H2 density^ in the range 2500 to 4000 cm~^. Given these 
ranges of temperature and density, we find that the ortho-H20/^^CO abundance ratio is a 
factor of 4 - 27 times as large as the H2O lio — loi / ^^CO J = 1 — hue ratio. The typical 
observed H2O lio — loi / ^^CO J = 1 — line ratio of 0.065 therefore implies an ortho- 
H20/^^C0 abundance ratio of 0.26 to 1.8. Assuming a ^^CO abundance of 10~^ relative 
to H2, a ^^CO/^'^CO abundance ratio of 30 for Galactic Center material (Langer & Penzias 
1990), and an ortho-to-para ratio of 3 for water vapor, we obtain a value of 1.2 - 8 xlO~^ 
for the typical total water abundance relative to H2 within the "180-pc" molecular ring. 

These water abundances inferred for the low density gas responsible for the Wlsr ~ 
80 — 120 km s^^ extended water emission are two to four orders of magnitude larger than 
those measured within cold dense molecular clouds (e.g. Snell et al. 2000). Thus the results 
obtained here from observations of water emission confirm our previous conclusion (NOO) - 
based upon absorption line observations - that the water abundance is typically larger in 
low-density interstellar gas. The most plausible explanation lies in the density dependence 



^For both molecules, rate coefficients are separately available for collisions with ortho- and para-H2. In 
the case of CO, the calculated rate coefficients are very similar for collisions with ortho- and para-H2, and 

thus the results of the statistical equilibrium calculation show little sensitivity to the assumed ortho-to-ratio 
ratio (OPR) for H2. For H2O, by contrast, the rate coefficients for excitation by ortho-H2 are typically 
an order of magnitude larger than those for excitation by para-H2; the derived water abundances therefore 
depend upon the assumed H2 OPR. Following Snell et al. (2000), we adopted the H2 OPR expected in local 
thermodynamic equilibrium at the gas temperature. 

*i.e. we assume that each coUisional excitation of the H2O lio — loi transition is followed by the emission 

of a line photon, even though the line is optically thick. 

^Towards 4 positions [(-3.2',-3.2'), (0.0',-9.6'), (0.0' ,-6.4'), and (0.0',- 3.2')], SWAS clearly detects ^^CO 
J = 5 — 4 emission in the 80 - 120 km s~^ range. Assuming a gas temperature of ~ 35 K and a plausible 
range of ^^CO column density (10^^ - 10^'^ cm'^), we find that the ^^CO J = 5 - 4/J =1-0 ratio is 
consistent with a H2 density of 2000 - 4000 cm~^. This is in good agreement with the range estimated by 
Kim et al (2002). 
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of the water vapor abundance, which is illustrated in Figure 6. This figure shows model 
predictions from the gas-grain model presented by Bergin et al. (2000) for three different 
molecular hydrogen densities. For lower densities the longer depletion timescales allow water 
vapor to exist in higher abundance. In dense regions, the abundance of water vapor is 
suppressed more rapidly by the depletion of oxygen nuclei onto icy grain mantles (Bergin et 
al. 2000). 

4.3. Water abundance in Sagittarius B2 

In a previous paper (NOO), we published the spectra of H2O and H2*0 Iqi — loi transitions 
obtained toward the (0, 0) position and shown here in Figure 7. Now that the general 
properties (i.e. velocity ccntroid and line width) of the extended 80 - 120 kms^^ emission 
have been elucidated, we return to a consideration of the H2O and H2^0 Iqi — loi absorption 
features at vlsr ~ 65kms~^. 

We have attempted to fit the SWAS-observcd absorption at the (0, 0) position with the 
use of two models. In the first model, we follow Zmuidzinas et al. (1995, hereafter Z95; 
see also Neufeld et al. 1997) in assuming that the absorbing water vapor is mixed with the 
continuum-emitting dust in the core of Sgr B2 itself; whereas in the second model, we assume 
the water vapor to he in separate intervening clouds - located close to Sgr B2 - that partially 
cover the continuum emission (see Ceccarelli et al. 2002). 

4.3.1. Water vapor in the core of Sgr B2 

In modeling the core of Sgr B2, we made use of a radiative transfer code very kindly 
made available by J. Zmuidzinas. This code yields a solution to the equations of statistical 
equilibrium for the water level populations - as a function of position within a spherically- 
symmetric, dusty, molecular cloud core - using the method of Accelerated Lambda Iteration 
(ALI) to obtain a converged, self-consistent solution to the level populations and the radiation 
field. Ray tracing is then used to determine the spectrum observed by a telescope with a 
Gaussian beam profile of specified projected diameter. The density in the core of Sgr B2 
was assumed to follow the profile of LG89, Model C, and the temperature of the gas and 
dust to follow the profiles assumed by Z95. The total continuum opacity was adjusted to 
match the observed (optically-thin) continuum emission at 550 GHz, with r assumed to 
vary with wavelength as A~^-^. We adopted the values given by PhiUips, Maluendes & Green 
(1996) for the rate coefficients for excitation of H2O in collisions with H2, assuming an local 
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thermodynamic equilibrium (LTE) ortho-to-para ratio for H2. The excitation model included 
the lowest 8 states of ortho-H20. 

In fitting the SWAS-observed H2^0 hne, we varied three parameters: the microturbulent 
hne width in the cloud core, the H2^0 abundance, and the velocity centroid relative to the 
LSR. All three were assumed to be constant throughout the cloud core. The blue curve 
in Figure 7 (lower panel) shows the best fit to the H2^0 absorption at ?;lsr ^ 65kms~^ 
obtained by this method. The best-fit microturbulent line width is 20kms^^, the ortho- 
H2^0 abundance is 6.4 x 10^^° relative to H2, and the LSR velocity centroid is 64kms~^. 
Assuming an H2^0/H2^0 abundance ratio of 250 (NOO; Langer & Penzias 1990), and a water 
ortho-to-para ratio of 3, we derive a total H2^0 abundance of 2.1 x 10^^ relative to H2. This 
derived abundance is in acceptable agreement with the value of 3.3 x 10~^ derived previously 
(Z95) from KAO observations of H^^O and reported^° by NOO. 

In the upper panel of Figure 7, the blue curve shows the fit to the H2^0 absorption at 
^LSR ~ 65kms~^, for the water abundance, line width and centroid velocity given above. 
One additional free parameter is introduced here: the peak antenna temperature (taken as 
0.18 K in this fit) of an emission component with the line centroid and line width derived 
in §3.4 above for the "180-pc molecular ring". The fit to the Hg^O absorption (upper panel. 
Figure 7, blue curve) is clearly deficient in three regards: (1) the predicted absorption line 
is too broad; (2) the predicted absorption line is too deep, and (3) the predicted emission 
peaks at vlsr = 33 and 95kms~^ are not observed. 

The first of these deficiencies is a direct consequence of the simplified velocity struc- 
ture assumed in this model. Observations of ammonia absorption towards Sgr B2 (e.g. 
Hiittemeister et al. 1995, hereafter H95) have revealed the presence of three distinct absorp- 
tion components in the t^LSR = 60 to 70 km s~^ interval (see discussion in the next subsection). 
Although none of these components shows a line width greater than ~ 12kms~^ (FWHM), 
collectively they blend to impart an overall line width ~ 20 km s^^ (FWHM) to optically-thin 
absorption lines. For the optically-thick Hg^O line, however, the observed line width does 
not increase as rapidly as it would for a single Gaussian component. 

Taken at face value, the second deficiency - the fact that the observed H2O line is less 



^°In order to fit the strengths of other, higher-excitation, H2*0 hnes detected toward this source, Neufeld 
et al. (1997) assumed that a higher water abundance of 5 x 10~^ was present in the inner parts of the Sgr B2 
core, whore the dust temperature was sufficient ( ^ 90 K) to cause the vaporization of water ice from grain 
mantles. Because the absorption in the lio — loi is completely dominated by material cooler than 90 K, the 
assumed presence or absence of enhanced water abundances in the region warmer than 90 K has no effect 
upon the fitting procedure adopted here. 
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deep than the model predicts - would imply that the Hg^O abundance is much less than 
assumed in the model. In Figure 8, middle panel, we show the line-center flux-to-continuum- 
fiux ratio - as predicted by the Zmuidzinas radiative transfer code - as a function of the 
assumed water abundance. The observed ratio ~ 0.25 would imply an Hg^O abundance of 
only ~ 2 X 10^^ and a very small H2^0/H2^0 ratio of ~ 3. We are aware of no chemical 
fractionation mechanism capable of yielding such a small H2^0/H2^0 ratio. A far more likely 
explanation of the anomalous depth of the H2^0 absorption feature is that some fraction 
(~ 25%) of the continuum emission detected in the SWAS beam originates in material with 
a very small water abundance (and which is not covered by any intervening cloud containing 
water). The likely presence of continuum emission that is not covered by a significant 
water column density has implications for the H2^0 abundance obtained above. For the 
water-covered component, the effective line-center to continuum ratio decreases from 0.49 to 
(0.49 — 0.25)/0.75 = 0.32, and the implied ortho-H2^0 abundance increases by a factor of 
almost 2 from 6.4 x 10"^'^ to 1.1 x 10"^. Given the same ortho-to-para and H2^0/H2^0 ratios 
assumed previously, this corresponds to a total H2^0 abundance of 3.7 x 10~^. Therefore, 
we conservatively adopt the range 2 — 4 x 10~^ as the H2^0 abundance implied by our model 
of the Sgr B2 core. 

The third deficiency of the simple model we considered is its prediction of narrow emis- 
sion peaks at vlsr — 35 and 95kms~^; this again is inconsistent with observed spectrum. 
A "double-horned" structure is predicted because the lio — loi resonance-fine photons that 
are excited coUisionally within the cloud core are scattered significantly in frequency before 
they can escape the cloud. The top and bottom panels in Figure 8 show the predicted fine- 
to-continuum ratios in these peaks and the predicted frequency offsets from line center (in 
velocity units). 

The observed absence of the blue-shifted horn is naturally explained by the known pres- 
ence of foreground absorbing material at i'LSR ~ 33kms~^. The absence of a red-shifted horn 
could easily be explained if the velocity dispersion increases slightly with increasing distance 
from the core center; such an increase would lead to the red-horn photons being spread over 
an spatially-extended region, leading to a decrease in the observed surface brightness relative 
to what is predicted by our constant velocity-dispersion model (see discussion in NOO). 



4.3.2. Water vapor in an intervening cloud 

As an alternative means of interpreting the H2O and Hg^O absorption towards the (0, 0) 
position in Sgr B2, we have fitted the spectrum with a series of absorption features that are 
assumed to originate in foreground molecular gas outside the core of Sgr B2. In this picture 
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- favored by H95 and by Ceccarelli et al. (2002; hereafter C02) for the origin of the NH3 

absorption observed toward Sgr B2 - the absorption originates in warm, low-density gas 
located outside the central core of Sgr B2 itself. According to C02, ISO observations of 
far-infrared NH3 absorption lines imply a gas temperature ~ 700 K, while the absence of 
observable far-IR CO emission places an upper limit n(H2) ,^ 3 x 10^ cm~^ on the gas density; 
in the best- fit model of C02, the warm absorbing gas is located at a distance ~ 1.15 pc from 
the center of the Sgr B2 core. 

In fitting the SWAS-obscrvcd H2O and Hg^O spectra in this picture, we adopted the 
absorption line parameters derived by H95 from high resolution observations of NH3 absorp- 
tion lines. Those observations were carried out separately toward two condensations within 
the Sgr B2 cloud: Sgr B2(N) and Sgr B2(M). These condensations, which are separated by 
~ 0.8', are not resolved individually by the 3.3' x 4.5' SWAS beam. H95 observed absorption 
components at i'lsr = 60.4 and 69.0 km s~^ toward Sgr B2(M) but not Sgr B2(N), and ab- 
sorption components at vlsr = 66.7 and 81.5 km s~^ toward Sgr B2(N) but not Sgr B2(M). 
The line widths and NH3 column densities inferred by H95 for these absorption components 
are listed in Table 3. 

We regarded the velocity centroids and fine widths for these four components to be 
fixed in our fitting of the SWAS-obscrved H2O and H2^0 spectra. We took as adjustable 
parameters the following: (1) the fraction, of the continuum emission covered by the 60.4 
and 69.0 km s~^ clouds^'^; (2) the fraction fs, of the continuum emission covered by the 66.7 
and 81.5 kms^^ clouds; (3) the strength of the 80 - 120 kms^^ emission feature at the (0,0) 
position in Sgr B2 (with an assumed line centroid and width as derived from the average 
spectrum shown in Figure 5 and discussed in §4.2); (4) the H2^0 line center optical depths 
in each of the four absorption components. Following C02, the absorption was considered to 
occur in a "foreground screen" . 

As in §4.3.1, we assumed an H2^0/H2^0 abundance ratio of 250, and a water ortho- 
to-para ratio of 3. In the upper panel of Figure 7, the red curve shows our best fit to the 
H2^0 absorption using this model. Because the H2^0 lines have optical depths of several 
hundred and therefore he on the fiattest part of the curve-of-growth, the Hl^O spectrum is 
almost entirely insensitive to the assumed optical depth; it therefore provides a relatively 
unambiguous determination of the covering factors J'a and Jb- The red curve represents our 
best fit to the H2^0 spectrum, obtained for = 0.6, /b — 0.2, and with a 80 - 120kms~^ 



^^Thc 60.4 and 69.0 km s ^ clouds were paired in this manner - with a common assumed covering factor 
- and the 66.7 and 81.5 km s~^ clouds likewise paired, because the H95 observations suggest that the former 
pair cover Sgr B2(M) but not Sgr B2(N), and the latter pair cover Sgr B2(N) but not Sgr B2(M) 
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emission feature of peak antenna temperature 0.2 K. 

Once the covering factors have been determined from the H2^0 spectrum, the hne-center 
optical depths can be determined from the H2^0 spectrum. The red curve in the lower panel 
of Figure 7 represents our best fit to the data, obtained with hne-center optical depth of 1.3, 
1.0, 0.8, and 1.5 respectively for the i'lsr = 60.4, 66.7, 69.0, and 81.5 kms~^ features. 

ISO observations of the 22i — lio transition of H2^0 at 108.073 /xm yield a further 
observational constraint upon the material responsible for the water absorption. A Long- 
Wavelength Spectrometer Fabry-Perot spectrum of the source (Goicoechea & Cernicharo 
2003) yields an absorption line equivalent width ~ 20kms~''^. Using a simple excitation 
model to estimate the fractional populations of the H2^0 rotational states, we find the 
observed 108 fim equivalent width to be consistent with the best-fit optical depths and 
covering factors derived from the SWAS observations (Table 3). Unfortunately, the 108 fim 
line lies on the flat part of the curvc-of-growth, so the observed equivalent width does not 
place strong constraints upon the physical conditions in the absorbing material. 

The absorbing column densities implied by the H2^0 line-center depths are inversely 
proportional to the difference /oi — /lo, where /oi is the fraction of ortho-water molecules in 
the loi rotational state and /lo is the fraction in Iiq. Thus the minimum column density of 
Hg^O is obtained for the assumption /oi = 1, /lo = 0; these minimum values are tabulated in 
Table 3. The beam-averaged, minimum total column density - weighted by the continuum 
antenna temperature - is A'"(H2^0) = 1.0 x 10^^cm~^, a value in excellent agreement with 
the value of 1.1 x lO^^cm"^ derived previously for the minimum Hg^O column density from 
KAO observations by Z95. The actual column density is larger by a factor (/oi — /io)~^- 

We have used a statistical equihbrium model for the Hg^O level populations to estimate 
the correction factor (/oi — /io)~^- In Figure 9, we plot the values of /oi, /lo and their 
difference /oi — /lo as a function of the radiation field. Here we follow C02 in assuming 
that the far-infrared continuum radiation originates in a region of radius -Rfir = 0.8 pc, 
and determine the H2^0 level populations as a function of the "dilution factor", W. The 
spectral shape of the far-infrared continuum source was obtained from Pipeline- 10 archived 
ISO LWS grating spectra and re-processed using the ISO Spectroscopic Analysis Package 
version 2.1.^^ The dilution factor is defined such that the sohd angle subtended by the 
far-infrared continuum source is AttW at the location of the absorbing water vapor. It is 
related to the distance, R, from the center of Sgr B2 (top horizontal axis in Figure 9) by 
the expression W — |[1 — (1 — Rfir/ R^y^^]- Solid curves in Figure 9 show results for 



^^The ISO Spectroscopic Analysis Package is a joint development by the LWS and SWS Instrument Teams 
and Data Centers. Contributing institutes are CESR, IAS, IPAC, MPE, RAL and SRON. 



-13- 



the temperature and density favored by C02 for the absorbing material: an H2 density of 
3 X 10^ cm~^ and a gas temperature of 700 K. Dotted curves show the results obtained for 
the same assumed temperature but with the density increased by a factor of 10. 

Two features of Figure 9 are noteworthy. First, the level populations are not greatly 
affected by a factor 10 increase in the assumed density. This behavior is explained by the fact 
that the gas density is very much smaller than the "critical density" for the Iiq— Iqi transition; 
the excitation of this transition is therefore controlled by the radiation field. Second, the 
quantity /oi — /lo is close to unity unless the absorbing material is located very close to 
the center of the Sgr B2 core. For the distance of 1.15 pc favored by C02, corresponding to 
W = 0.15, the value of foi — fio is 0.8 and the correction factor (/oi — /lo)^^ is only 1.25. 
Significantly smaller values of the distance of the absorbing material from the cloud core are 
ruled out by the requirement that the absorbing water vapor must lie in front of most of 
the emitting dust. We conservatively adopt an upper limit of 1.5 for the correction factor 
(/oi - /io)~^- 

Using the analysis method of §4.3.2, where the water absorption is assumed to arise in 
a foreground cloud of warm, low-density material, our estimate of the beam-averaged H^^O 
column density lies in the range 1.0 — 1.5 x 10^^ cm~^. Given an assumed abundance ratio 
of 250 for H2^0/H2^0, the Hg^O column density we derived implies an H2^0 column density 
of 2.5 — 4 X 10^^ cm^^. The total H2 column density of the warm absorbing gas layer is not 
well constrained by observations; if we adopt the estimate of 10^^ cm~^ given by Flower et 
al. (1995), then the water abundance is a few xlO~^. 

As noted by Flower et al. (1995) - and contrary to the assertion of C02 that there is 
no discrepancy - the observed water column density lies an order of magnitude below the 
predictions of shock models for the warm, absorbing gas layer in front of Sgr B2. Recent ob- 
servations of high OH abundances in Sgr B2 (Goicoechea & Cernicharo 2002) may provide an 
explanation of the discrepancy between theory and observation. As suggested by Goicoechea 
& Cernicharo, if the shocked gas is irradiated by strong ultraviolet radiation from a popu- 
lation of hot stars, then photodissociation will reduce the H2O abundance and increase the 
OH abundance. Further chemical modeling is needed, however, to determine whether this 
explanation is either quantitatively plausible or consistent with the large observed abundance 
of NH3 in the absorbing gas. 



^"^Wc note that this vahic is substantiaUy smaUcr than an estimate of 10^^ cm ^ given by C02 for A''(H2*0) 
but agrees with that obtained by Comito et al. (2003). 
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4-3.3. Summary: water vapor in the cold core or in a warm intervening cloud? 

In §4.3.1 and §4.3.2 above, we have used two alternative methods to analyse the SWAS 
observations of H2^0 and Hg^O at the (0, 0) position in Sgr B2. Although the first method 
- in which the water vapor and dust are assumed to be well mixed in the cloud core - yields 
a fit to the data which is inferior to that of the second method - in which the water vapor 
is assumed to he in a foreground gas cloud - the poorer fit can be ascribed to the fact that 
the assumed velocity structure was far simpler in the first analysis method. On the basis 
of the SWAS data alone, therefore, we do not find compelling evidence to favor one picture 
over the other. 

In principle, both the cold core and a warm intervening cloud might contribute to the 
observed water absorption. Insofar as the water column densities involved add linearly, we 
may conclude that 

/ iVfg(H20) X / Xc(H20) X ^ ^ 

V2.5to4 X 10i6cm-2y V2-4xl0-7; ^' 
where A^fg(H20) is the water column density in the warm foreground material and a;c(H20) 
is the water abundance relative to H2 in the core. Clearly, we may regard 4x10^^ cm"^ and 
4 X 10~^ as conservative upper hmits on A^fg(H20) and Xc(H20).^^ 
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-'^''We note, however, that Comito et al. (2003) have carried out an analysis which combines observational 
data from SWAS with ground-based observations of HDO and para-Hj^O. They conclude that a significant 
fraction of the water absorption must take place within warm foreground material. 

Similar remarks apply to the far-IR hydrogen fluoride absorption observed by Neufeld et al. (1997) 
toward Sgr B2. Assuming the HF absorption to occur entirely in the core of Sgr B2, Neufeld et al. derived a 

low HF abimdance that implied a substantial depletion factor ^ 50 for fluorine in the dense core. C02 pointed 
out (correctly) that a considerably larger HF abundance would be implied if the absorption were assumed 
to take place in intervening material at a distance ~ 1.15 pc from the core center, but argued (erroneously) 
that this would weaken the argument for a high fluorine depletion in the core of Sgr B2. Clearly, if part 
or all of the observed HF absorption occurs in foreground material outside the core of Sgr B2, then the HF 
abundance within the core would have to be even less than that inferred by Neufeld et al. (1997). This 
would merely strengthen the argument for a high fluorine depletion in the core of Sgr B2. 
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Fig. 1. — SWAS spectra of the 556.936 GHz lio — loi transition of H2^0, observed in the vicinity of 
Sgr B2. Each panel covers the LSR velocity range —145 to 145 km s^^. Adjacent spectra are offset 
by 3.2'. Horizontal dotted lines indicate the the continuum flux level and the zero flux level. The 
vertical dotted lines indicate a LSR velocity of 65kms~^. On each panel, a pair of numbers, X/Y, 
describe, in order, the continuum antenna temperature and the antenna temperature (in K) at the 
top of the panel. The feature at ^lsr ~ ~90 to — 50kms~^ is a ^^CO emission feature detected in 
the lower sideband. The (0,0) position is in the middle of the central group of nine spectra. 
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Fig. 2. — SWAS spectra of the 550.927 GHz J = 5 - 4 transition of ^^CO, observed in the vicinity 
of Sgr B2. Each panel covers the LSR velocity range —145 to 145kms~^. The vertical dotted lines 
indicate a LSR velocity of 65kms^^. Adjacent spectra are offset by 3.2'. The number on each 
panel gives the antenna temperature (in K) at the top of the panel. The feature at vlsr ~ —120 
to -80kms-i is H2O emission being detected in the upper sideband. The (0,0) position is in the 
middle of the central group of nine spectra. 
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Fig. 4.— SWAS RfO spectrum of Sgr A, with the beam centered at a = 17'' 45™ 39. '7, 
S = -29°00'30" (J2000), corresponding to an offset (-21.8', -37.4') relative to Sgr B2. This 
spectrum was obtained in 454 hours of total integration time at 11 separate local oscillator 
(LO) settings during the months of July 1999, April 2000, May 2000, August 2000, and Sept. 
2000. The '^^CO J=5-4 and H2^0 lio — loi lines were detected simultaneously in different 
sidebands of SWAS Receiver 2 (see text). The use of different LO frequencies allowed us to 
reconstruct the true ^^CO J=5-4 line profile free from contamination by water emission or 
absorption. This profile was then subtracted from the Hg^O spectrum (with proper account 
taken of the varying spacecraft velocity and different LO settings used), yielding a H2^0 
spectrum free of ^^CO emission which has a considerably larger bandwidth (~ 500 km s~^) 
than a typical SWAS spectrum (~ 210 km s~^). The rms noise in the spectrum varies from 
0.015 K for t>LSR in the range -60 to 300 km s~^ to 0.05 K for velocities outside that range. 
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Fig. 5. — Average H2^0 Iiq — loi spectrum for those 14 positions that are devoid of absorption 
at the Sgr B2 systemic velocity. The sohd hne is a Gaussian fit to the average emission 
feature with the following parameters: velocity centroid 93kms~^, line width 43kms~^ 
(FWHM), peak antenna temperature 0.075 K, and velocity-integrated antenna temperature 
3.4Kkms^^. The observed emission is believed to orignate in the far side of the "180-pc" 
molecular ring. 
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Fig. 6. — Water vapor abundance as a function of time for three different H2 densities. Model 
results are predicted using the gas-grain model with limited surface chemistry presented by 
Bergin et al. (2000). The gas is assumed to be completely shielded from external ultraviolet 
radiation and to have a temperature of 10 K. 
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Fig. 8. — Predictions of the Zmuidzinas radiative transfer code for Sgr B2. Top panel: peak 
line-to-continuum ratio for the double horns, as a function of the total water abundance. 
Middle panel: ratio of line-center flux to continuum flux, as a function of the total water 
abundance. Lower panel: velocity separation ofthe double horns from line center. 
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Fig. 9. — Fractional level populations, /oi and /lo, and their difference /oi — /lo, as a function 
of distance from the center of Sgr B2 (top horizontal axis; the bottom horizontal axis shows 
the corresponding "dilution factor" for the radiation field [see text]). Solid curves apply 
to an H2 density of 3 x 10^ cm~^ and a gas temperature of 700 K, values favored by C02. 
Dotted curves show the results obtained for the same assumed temperature but with the 
density increased by a factor 10; the results are virtually unchanged because the molecular 
excitation is controlled by the radiation field, not coUisional excitation. 
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TABLE 1 

H2^0 observations toward Sgr B2: inferred continuum fluxes 



Offsef* Time ^ Continuum T4 550 GHz continuum 350 /xm continuum 
(min) (mK) flux'^ (kJy) flux*^ (kJy) 



(-9.6, +0.0) 


394 


50 


0.52 


0.06 


(-6.4, +0.0) 


332 


100 


1.04 


1.71 


(-3.2,-3.2) 


157 


161 


1.67 


7.49 


(-3.2, +0.0) 


107 


292 


3.03 


10.8 


(-3.2, +3.2) 


263 


181 


1.88 


4.90 


(+0.0,-9.6) 


404 


< 20 


< 0.3 


1.77 


(+0.0,-6.4) 


495 


64 


0.66 


3.27 


(+0.0,-3.2) 


99 


332 


3.45 


12.7 


(+0.0, +0.0) 


373 


618 


6.42 


31.9 


(+0.0, +3.2) 


100 


221 


2.29 


11.0 


(+0.0, +6.4) 


487 


20 


0.21 


1.02 


(+0.0, +9.6) 


391 


45 


0.47 


0.03 


(+3.2,-3.2) 


263 


107 


1.11 


3.82 


(+3.2, +0.0) 


102 


337 


3.50 


8.32 


(+3.2, +3.2) 


261 


95 


0.99 


5.68 


(+6.4, +0.0) 


202 


105 


1.09 


1.78 


(+9.6, +0.0) 


379 


94 


0.98 


0.10 


(+12.8, +0.0) 


394 


95 


0.99 


<0.025 


(+16.0, +0.0) 


225 


80 


0.83 


<0.025 



" Offset {a cos 6, 6) in arcmin relative to a = 17^ 47™ 19.^7, S = -28°23'07" (J2000) 
^ On-source integration time 

Antenna temperature at 550 GHz 

Measured by SWAS 
^ Measured by SHARC (Dowell et al. 1999) 
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TABLE 2 

Hg^O observations toward Sgr B2: emission line strengths 



Offset" Jso {Tl - TXc)dv HaO/i^CO 
(K kms~^) line ratio*^ 



(-9.6, +0.0) 


0.74 


0.074 


(-6.4, +0.0) 


1.07 


0.046 


(-3.2,-3.2)'^ 


1.66 


0.035 


(-3.2, +0.0)'^ 


3.33 


0.053 


(-3.2, +3.2) 


1.50 


0.046 


(+0.0,-9.6) 


< 0.76^ 


< 0.03P 


(+0.0,-6.4) 


2.15 


0.043 


(+0.0, -3.2)^^ 


2.89 


0.044 


(+0.0, +0.0)<^ 


4.18 


0.103 


(+0.0, +3.2)<^ 


2.31 


0.071 


(+0.0, +6.4) 


3.04 


0.060 


(+0.0, +9.6) 


1.98 


0.058 


(+3.2,-3.2) 


1.49 


0.036 


(+3.2, +0.0) 


3.50 


0.070 


(+3.2, +3.2) 


3.53 


0.071 


(+6.4, +0.0) 


4.09 


0.098 


(+9.6, +0.0) 


3.34 


0.125 


(+12.8, +0.0) 


2.43 


0.101 


(+16.0, +0.0) 


2.16 


0.089 



" Offset {a cos 6,6) in arcmin relative to a = 17^47"^ 19.^7, 6 = -28°23'07" (J2000) 
Continuum-subtracted integrated over the interval ^lsr = 80 to 120 kms~^ 
Ratio of continuum-subtracted T|, integrated over the interval wlsr = 80 to 120 

kms~^, for the H2O lio — loi transition and the ^^CO J = 1 — transition 

(Bally et al. 1987). This ratio is computed using fluxes that are uncorrected for 

antenna efficiency because the efficiencies of both antennas are comparable 

(Bally, priv. comm.; Melnick et al 2000). 
Line ratio computed for the smaller interval vlsr = 93 to 120 kms~^ so as to 

avoid the influence of water absorption near the Sgr B2 systemic velocity (see text). 

^ 3 sigma upper limit 
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TABLE 3 
Water absorption line parameters 



LSR velocity'* Line width" Covering t(H2^0) Minimum 
(kms-^) (kms-^FWHM) factor A^(H^^O) 



60.4 8.0 0.6 1.3 1.7 xl0i^cm-2 
66.7 12.1 0.2 1.0 2.0 X 10^6 cm-2 
69.0 8.0 0.6 0.8 LlxlO^^cm'^ 

81.5 9.8 0.2 1.5 2.5 xl0^^cm-2 



" From H95 

^ Hl^O optical depth at line center 



